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Abstract 

We elaborate on a recently proposed extension of the Gerasimov-Drell-Hearn (GDH) sum rule 
which is achieved by taking derivatives with respect to the anomalous magnetic moment. The new 
sum rule features a linear relation between the anomalous magnetic moment and the dispersion 
integral over a cross-section quantity. We find some analogy of the linearized form of the GDH 
sum rule with the 'sideways dispersion relations'. As an example, we apply the linear sum rule to 
reproduce the famous Schwinger's correction to the magnetic moment in QED from a tree-level 
cross-section calculation and outline the procedure for computing the two-loop correction from a 
one-loop cross-section calculation. The polarizabilities of the electron in QED are considered as well 
by using the other forward-Compton-scattering sum rules. We also employ the sum rules to study 
the magnetic moment and polarizabilities of the nucleon in a relativistic chiral EFT framework. In 
particular we investigate the chiral extrapolation of these quantities. 
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I. INTRODUCTION 



The Gerasimov, Drell, Hearn (GDH) sum rule (SR) [1], which relates a system's anoma- 
lous magnetic moment to a weighted integral over a combination of doubly polarized pho- 
toabsorption cross sections, has received a good deal of attention in recent years. Impressive 
experimental programs to measure these photoabsorption cross-sections for the nucleon have 
recently been carried out at ELSA and MAMI (for a review see Ref. [2]). Such measurements 
provide an empirical test of the GDH SR, and can be used to generate phenomenological 
estimates of electromagnetic polarizabilities via related SRs, as will be discussed below. The 
GDH SR is particularly interesting because both its left- and right-hand-side can be reli- 
ably determined, thus providing a useful verification of the fundamental principles (such as 
unitarity and analyticity) which go into its derivation. At the present time the proton sum 
rule is satisfied within the experimental precision, while the case is still out for the neutron. 
However, it is not the purpose of this paper to discuss these experimental aspects of the SR, 
but rather to see what can be learned on the theoretical side. 

Recently, we have shown [3] that by taking derivatives of the GDH sum rule with respect 
to the anomalous magnetic moment one can obtain a new set of sum rule-like relations 
with intriguing properties. In particular, this procedure provided a sum rule involving 
the anomalous magnetic moment linearly rather than quadratically, which allows for the 
derivation of quantities such as the Schwinger moment in a much simpler fashion than via 
the usual GDH method. In this paper, we shall further examine these forms and apply them 
to the nucleon magnetic moment. 

After a lightning review of the GDH and related sum rules (Sect. II), we derive the 
modified versions (Sect. HI) and demonstrate that some of these relations have the form 
of so-called sideways dispersion relations (Sect. IV). Then we consider applications to the 
nucleon in the context of chiral perturbation theory and show how the new sum rules allow an 
elementary calculation (to one loop) of quantities such as magnetic moments (Sect. V) and 
polarizabilities (Sect. VI) to all orders in the heavy baryon expansion. The chiral behavior 
of the nucleon magnetic moments and polarizabilities is addressed in Sect. VII. Returning to 
QED in Sect. VIII, we demonstrate how the new sum rules can be applied to the two- loop 
calculation of the anomalous magnetic moment in a straightforward fashion. In the final 
section, we summarize our findings and suggest prospects for future work. 

II. COMPTON-SCATTERING SUM RULES IN QED 

The forward-scattering amplitude describing the elastic scattering of a photon on a target 
with spin s (real Compton scattering) is characterized by 2s+l scalar functions which depend 
on a single kinematic variable, e.g., the photon energy cu. In the low-energy limit each of 
these functions corresponds to an electromagnetic moment — charge, magnetic dipole, electric 
quadrupole, etc. — of the target. In the case of a spin-1/2 target, the forward Compton 
amplitude is generally written as 

T{uj) = t* ■ eficu) + ia- {t* x e) g{uj) , (1) 

where e' is the polarization vector of the incident and scattered photon, respectively, 
while (T are the Pauli matrices representing the dependence on the target spin. The crossing 
symmetry of the Compton amplitude of Eq. (1) means invariance under e' e, cu —cu, 
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which obviously leads to f{uj) being an even and g{uj) being an odd function of the energy — 
f{uj) = f{—uj), g{ijj) = —g{—u)). The two scalar functions /(<^), g{oj) admit the following 
low-energy expansion, 



f{^) = + {o^E + /3m) uj' + , (2a) 

^^^^ ^ -^^ + ^°^' + ^^^')' (2b) 

and hence, in the low-energy limit, are given in terms of the target's charge e and anomalous 
magnetic moment (a.m.m.) k. The next-to-leading order terms are given in terms of the 
nucleon electric (a;^;), magnetic (Pm), and forward spin (70) polarizabilities. 

In order to derive sum rules (SRs) for these quantities one assumes the scattering ampli- 
tude is an analytic function of u everywhere but the real axis ^, which allows writing the real 
parts of the functions /(a;) and g{u!) as a dispersion integral involving their corresponding 
imaginary parts. The latter, on the other hand, can be related to combinations of doubly 
polarized photoabsorption cross-sections via the optical theorem, 

Im /(a;) = ^ [(T1/2H + 0-3/2(0;)] , (3a) 

Im g{uj) = ^ [(7i/2(w) - 0-3/2 (w)] , (3b) 

where is the doubly-polarized total cross-section of the photoabsorption processes, with 
A specifying the total helicity of the initial system. Averaging over the polarization of initial 
particles gives the total unpolarized cross-section, ctt = \{p\ii + c-3/2). 

After these steps one arrives at the results (see, e.g., [4] for more details): 

2 pOO ( f\ 







with Act = (T3/2 — cri/2, and where the sum rule for the unpolarized forward amplitude /(<^) 
has been once-subtracted to guarantee convergence. These relations can then be expanded 
in energy to obtain the SRs for the different static properties introduced in Eq. (2). In this 
way we obtain the Baldin SR [5, 6]: 

1 ariuj) , 
o:e + Pm = -z-^ I 5— duj , (6) 



the GDH SR: 



a SR for the forward spin polarizability: 



1 r A^(^) , 



Resonance poles may occur but lie on the second Riemann sheet. 
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and, in principle, one could continue in order to isolate higher order moments [7]. 

Let us now see how these sum rules are of use in field theory. We first consider the 



case of the electron in QED. To lowest order in the fine-structure constant, ctg 



V47r, 



the photoabsorption process is given by the tree-level Compton scattering. The tree-level 
helicity amplitudes for this process are well known [8] 
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where M is the electron mass; s, t and u arc the Mandelstam variables. Using these results 
the unpolarized and double-polarized cross sections can be determined as: 
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where we have defined x — uj/M. Substituting the latter expression into the r.h.s. of the 
GDH SR, Eq. (7), one can easily see that the integral vanishes exactly, which is required 
because otherwise the sum rule would lead to a nonsense result — the electron a.m.m. would 
receive contributions of order ^/owi". 

At next order, 0(0;^^), the l.h.s. of the GDH SR receives a nonzero contribution in the 
form of the Schwinger correction: k — aem/STT. In order to check that the same result is 
obtained on the r.h.s. of the sum rule is quite a formidable task since at this order one 
must know the Compton scattering amplitude to one loop, as well as account for the pair- 
production channel. Nevertheless, the calculation of the relevant helicity amplitudes was 
carried out more than three decades ago by Milton, Tsai, and deRaad and application to 
the GDH SR has relatively recently been performed by Dicus and Vega [9], who verified 
(numerically) that the GDH SR holds at O(aL) in QED- 

Recently we have found a much simpler method by which to verify the GDH SR at this 
order. This technique is briefly described in the next section, while in the rest of the present 
section we examine the sum rules for polarizabilities, Eqs. (6) and (8). At first glance, 
such sum rules do not appear to be of much utility, since when evaluated in the case of 
QED, the static polarizabilities obtained therefrom diverge. However, this is not a problem, 
as pointed out by Llanta and Tarrach [10], who emphasized that one can determine well 
defined asymptotic forms in the limit as uj approaches zero. Thus, for example we can 
perform the integration for the spin even/odd amplitudes and determine that 
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One can now use these expressions to define quasi-static polarizabilities, where by this we 
mean that these quantities contain both a constant term and a term which behaves as 
logo; — it is this latter piece that is responsible for the divergence as the static limit is taken. 
Thus, for example, we may define quasi-static values 
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for the sum of electric and magnetic polarizabilities and 
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(14) 



(15) 



for the forward spin polarizability. In this way we can also generate generalized sum rules 
for quadrupole and higher polarizabihties via Eqs. (4) and (5) [7]. 

It is interesting to note that it is possible to determine the nonanalytic component of these 
quasi-static moments in a simpler fashion — by using only the the low energy expansion of 
the cross sections. That is, while the integrals 

io r-a;2 

are A-dependent if n is an even integer, in the case that n is odd there exists a A-independent 
logarithm 
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which then determines the logarithmic component of the quasi-static polarizabilities. That 
is, since 
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we require that the nonanalytic piece of the polarizabilities have the form 
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in agreement with the exact forms given above. Likewise higher order forms can be deter- 
mined. 

We see then that the use of dispersive techniques in QED allows a straightforward extrac- 
tion of information about polarizabilities and about the anomalous magnetic moment. As 
expected, the latter is complete agreement with the result obtained by conventional means, 
while the former requires such dispersive methods, since the corresponding static quantities 
are divergent. 
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III. DERIVATIVES OF THE GDH SUM RULE 



We now review the derivation of the new form of sum rule. We begin by introducing a 
'classical' (or 'trial') value of the electron a.m.m., kq. At the Lagrangian level this amounts 
to the introduction of a Pauh term for the spin- 1/2 field : 

>Cpauli = ^V'^M.V'^"^ (21) 

where F^^'^ is the electromagnetic field tensor and a^^, = (i/2)[7^, 7,^] is the usual Dirac tensor 
operator. At the end of our calculation, we will set Kq to zero, but for now the total value 
of the a.m.m. is k, = + Sk, with Sk, denoting the loop contribution. (It is important to 
note that both 6k, and the cross-section become explicitly dependent on kq.) We then start 
taking derivatives of the GDH SR with respect to kq, which is subsequently set to zero, so 
that the total a.m.m. k returns to its usual loop value. We find 

Uir^aejM^) kk' = H ^a\uj) — , (22) 

Jo ^ 

{47:^aem/M^){K'^ + KK") = H Aa" (cu) — , (23) 

Jo ^ 

and so on. To lowest order in ctem we find 

{in'a,jM')nK^^-'^ = H Aa^-\u) — , (24) 

Jo ^ 

where k'^") denotes the nth derivative of k with respect to kq. This allows in principle 
the computation of k to order a"^ by using the 1st to nth derivatives of the cross-section 
computed to order ctg^^ to ctg^, respectively. 

In particular, to lowest order we have the result : 



The striking feature of this sum rule is the linear relation between the a.m.m. and the 
(derivative of the) photoabsorption cross section, in contrast to the GDH SR where k, appears 
quadratically, and although the cross-section quantity Act' (a;) is not an observable, it is very 
clear how it can be determined within a specific theory. Thus, for example, the first derivative 
of the tree- level cross-section with respect to Kq, at Kq — 0, in QED was worked out in [3]: 



9-7rrv^ 



2Muj / 3M\ , / 2uj 

6 - - 2 + ln l + — 

M + 2a; 2 \ ^ uj J \ M 



(26) 



It is not difficult to find then that 



oc 


Substituting this result in the linearized GDH SR, Eq. (25), we obtain k, = Q;em/27r — 
Schwinger's one-loop result. We emphasize that this result is reproduced here by computing 
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only a (derivative of the) tree-level Compton scattering cross-section and then performing 
an integration over energy. This is definitely much simpler than obtaining the Schwinger 
result from the GDH SR directly [9], which requires input at the one-loop level. Along these 
lines, however, one can facilitate the two- and more loop calculations. We elaborate on this 
possibility in Sect. VIII. 



IV. CONNECTION TO SIDEWAYS DISPERSION RELATIONS 

It is interesting to observe that, by changing the integration variable to s = -|- 2Muj, 
the linearized GDH SR, Eq. (25), can be written as 



K — 



47r2 



ds 



S-M2 



(28) 



M2 



which is a special case of a 'sideways dispersion relation' as is demonstrated below. 

The sideways dispersion for the a.m.m. is obtained by considering the half-off-shell elec- 
tromagnetic vertex: 



2M 
{p + p')^ 
2M 



(29) 



with q = p' — p the photon momentum, p' the off-shell particle momentum (s' = p'^), and p 
the on-shell momentum (p^ = M^). Furthermore, 



2M ' 



(30) 



are the positive- and negative-energy state projections, and F^"^^ and G^^^ are corresponding 
half-o£F-shcll form factors. In this decomposition of the vertex, the a.m.m. is identified as 

K = -G^+\M\0). 

Akin to the Compton amplitudes considered above, the form factors must be analytic 
functions of s' everywhere in the complex plane except along the unitarity cut, which lies 
on the real axis for Res' > M^. This analyticity allows a dispersive representation for each 
of the form factors: 



F^^\s',q') = IJ 
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and these relations are called sideways dispersion relations [11, 12]. 



7 



In particular, from the relation for we find 



oo 

1| ImG<^ (33) 

n J s- ' ^ ^ 

M2 

which is of the same form as the linearized GDH sum rule, Eq. (28). Specifically, one can 
show that 

ImGW(., 0) = -— Aa'{s)l^^, + — [Aa'{s) - Aa'Ml^^, , (34) 

where the second term integrates to 0. The function Acr^ is defined as the interference of the 
tree-level amplitude. Fig. 1, and the s-channel graphs with a Pauli-coupling insertion (first 
two graphs in Fig. 2). Recall that the function Act' is the interference of the tree-level graphs 
and both the s- and w-channel graphs with a Pauli-coupling insertion (the four graphs in 
Fig. 2), and analogous relations for the nucleon case will be discussed in the next section. 



V. SUM RULES IN CHIRAL PERTURBATION THEORY 

Consider now a theory of nucleons interacting with pions via pseudovector coupling: 

C^NN^^^ri'r'^^d.n'^, (35) 

where g is the pion- nucleon coupling constant, M is the nucleon mass, r" are isospin Pauli 
matrices, ip is the nucleon field and tt" is the isovector pion field. For our purposes this 
Lagrangian is sufficient to obtain the leading order results of chiral perturbation theory. 

To lowest order in the coupling g, the photoabsorption cross section in this theory is 
dominated by the single pion photoproduction graphs as displayed in Fig. 3(a), and we find 
the corresponding helicity difference cross sections : 

Aa^-"p\u;) = l{2as + l-x)\n^^-2X[x(a-2)+s{a + 2)]\, (36a) 

M'^x'^ I a — A ) 

^'"'"'("' = + (36b) 

Aff<"°"'(w) = 0, (36c) 

^<'"""('^> - |Vb.f±A + (2..-l-.)L.^-2.A}, (36d) 

which are expressed here in terms of the dimensionless quantities: 

C = {eg/A7rf , (37a) 

x^uj/M, iJL^m^/M, s^s/M^ = l + 2x, (37b) 

a^{s + M^ -ml)/2s, (37c) 

(3^{s-M^ + ml)/2s = 1 - a, (37d) 

A = (l/2s)y/s - (M + m^) -(M- m^f . (37e) 
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with rriTr denoting the pion mass. 

The nucleon anomalous magnetic moment is generated from loop diagrams and hence 
begins at 0{g'^), implying that the Ihs of the GDH SR begins at 0{g'^). Since the tree-level 
cross sections above are 0{g^), we must require, as in the case of QED, the consistency 
conditions 

oo 

^ A(t(^) (cu) = 0, for / = ttV, TT+n, 7r°n, Tf-p, (38) 



I. 



where a;th = rriT^il + mT,/2M) is the threshold of the pion photoproduction reaction. This 
requirement is indeed verified for the expressions given in Eq. (36) — the consistency of GDH 
SR is maintained in this theory for each of the pion production channels. 

We now turn our attention to the linearized GDH sum rule — Eq. (25). In this case we 
introduce Pauli moments and Kon for both the proton and neutron, respectively. The 
dependence of the resulting cross-sections on these quantities can then generally be presented 
as: 

Acr(a;; ^op, /ton) = Acr(a;) + KopAaip{uj) + Acrin(a;) + ... , (39) 
where we denote 



(40) 

K0p=K0n=0 



Furthermore, we introduce total proton and neutron photoproduction cross sections Aa*-^^ 
and Acr*^") and write the corresponding GDH SRs and their first derivatives. In this way we 
obtain the relations: 

(i) the GDH SRs: 



oo 



M2 '^P n .1 u 



(ii) the linearized SRs (valid to leading order in the coupling g) 



<^th 

oo 

= -J-Aa[J. (41c) 
(iii) the consistency conditions (valid to leading order in the coupling g): 

oo oo 

= i/^Ae 0=i/!^A.i:> (41d) 

71 J Ul 71 J UJ 
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The first derivatives of the cross-sections that enter in Eq. (41), to leading order in 
arise through the interference of Born graphs Fig. 3(a) with the graphs in Fig. 3(b) and we 
find: 

Aag) = Aaif + Aaff = {2x\[A + (1 - 2a) (1 + 2s)] - 2sA(a + 2) 

Ivl X 



- jji^x In ^ + (2as + 1 - a;) In \ , (42a) 

p — X a — X) 

AaiJ = AaL + Aa^ " m% |2A + /. In ^ - In --^ j , (42b) 
A.£).A.{r^) + A.{r"^ = |[^|ln^ + 2A(./3-.-a)}, (42c) 
Aa£) ^ Aai^°") + Aai;--) = 1^ I {2sa - x) In ^ + 2A(x - 2.") } . (42d) 

Using the latter two expressions we easily verify the consistency conditions given in Eq. (41d). 
Employing the linearized SRs, we obtain for the terms independent of n^p and «;on- 



oo 



' 1 - ^ (4- ll/^' + 3/ f^ ^^^^^g ^ _ 6;,2 ^ 2/^2 (-5 + 3 /x^) In /X i , (43a) 



(47r)2 1 yr^^ 2 



oo 



We have checked that the expressions of Eq. (43) agree with the one-loop calculation done 
by using the standard Feynman-parameter technique. To this order, the pseudoscalar pion- 
nucleon coupling gives exactly the same result, which can easily be verified by using the 
expressions of Appendix B. 

We note that the chiral expansion of the corrections Eq. (43) begins with a constant, 
i.e., at 0{1), and not with as is inferred by power counting. However this is not a 
problem, because the power-counting-violating constant term is absorbed into the counter 
term kq. The rest of the loop contribution obeys power counting. In general, in relativistic 
ChPT the power counting applies only to renormalized loop contributions, where the naively 
power-counting violating terms are absorbed into counter terms [13]. 

Even after such a renormalization of the constant terms in Eq. (43), our sum-rule result 
is no^ in agreement with the covariant ChPT calculation of Kubis and Meifiner [14], which is 
based upon the infrared regularization procedure of Becher and Leutwyler. The differences 
appear only in the terms that are analytic in the quark mass, rUq ~ m^. This discrepancy 
can be traced back to the fact that the 'infrared-regularized' loop amphtudes do not satisfy 
the usual dispersion relations — their analytic properties in the energy plane are complicated 
by an additional cut due to an explicit dependence on a/s. In other words, they do not obey 
the analyticity constraint which is imposed on the sum rule calculation. Certainly, as the 
differences are analytic in the quark mass, they can be reconciled order by order due to the 
inclusion of the appropriate counter terms. 
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The interpretation of the hnearized GDH sum rule is clarified by observing that it has the 
form of a sideways dispersion relation. As remarked above for the QED case, one can write 
down the following sideways dispersion relation for the loop contribution to the nucleon 
a.m.m.: 



oo 



, If. ImGW(s,0) 

= " J s-M^ • 

where is now the appropriate half-off- shell form factor for the nucleon, cf. Eq. (29). 
The absorptive part of this form factor to one loop is obtained by computing the indicated 
cuts of the graphs in Fig. 4. For instance, the result for the proton can be written as: 



P g2 



(45) 



where the first term is equivalent to the contribution of Fig. 4 (a), whereas the second term 
represents Fig. 4 (b). Explicit forms of these functions AcTg are listed in the Appendix B. 
Making the connection with the linearized GDH sum rule, we find that these AcXs can equiva- 
lently be determined by computing the interference of the Born pion-photoproduction graphs 
(Fig. 3) and the s-channel graph with a Pauli-coupling insertion [first graph in Fig. 3(b)]. 

Note that Acr^^p"'^ entering the above sideways dispersion relation is exactly the same 
as the corresponding integrand Acr^p entering the linearized GDH sum rule of Eq. (41b). 
The term Aa^JlJ'^ in the integrand of Eq. (44) on the other hand is different from the cor- 
responding term Acrjp in Eq. (41b). However, one can easily check that upon integration 
both integrands give the same contribution to the magnetic moments. We have verified this 
result for both pseudo-scalar and pseudo-vector ttNN couplings. 

In Appendix C, we, for completeness, list all the expressions for the second derivatives of 
the GDH cross-sections for different single-pion production channels in Born approximation. 

VI. CHIRAL CORRECTIONS TO NUCLEON POLARIZABILITIES 

It is well-known that the Heavy-Baryon ChPT (HBChPT) at order 0{p^) yields the 
prediction for the electric and magnetic polarizabilities of the nucleon: 

(HBLO) ^ STTOem ( _9A_\ ' ^ ^ ^ ^q_4 ^^3 ^^g^) 

^^HBLO) ^ ^ fgAV ^ ^ 2 X 10-^ fm^ , (46b) 

where qa — 1-26, /^r ~ 92.4 MeV. Here the couphngs are related to the ttNN couphng 
constant used in the previous section via the Goldberger-Treiman relation: qa/ fir = g/M. 

Eqs. 46 are a true prediction of 0{p^) HBChPT (there are no counter-terms at this order) 
and turn out to be in remarkable agreement with experiment, e.g., for the proton [15]: 



ij''^^ = 12.1 ± 0.3 (stat) ^ 0.4 (syst) ± 0.3 (mod) [x lO"^ fm^] , (47a) 
((^^p) ^ i Q ^ 0.4 (stat) ± 0.4 (syst) ± 0.4 (mod) [x 10"^ fm^] . (47b) 
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By using forward sum rules, we can discuss here only the sum of the electric and magnetic 
polarizabilities. On the experimental side, a recent determination from the Baldin's sum 
gives [16]: 

("E + PmYp"" = (13.69 ± 0.14) X 10-^ fm^ , 

(aE + PmYu^ = (14.40 ± 0.66) x 10"^ fm^ , (48) 

for proton and neutron, respectively. 

In order to find the leading order relativistic prediction of chiral loops we have com- 
puted the unpolarized total cross-sections, corresponding to the Born graphs of single-pion 
photoproduction : 



TlC 



[uj^ - f/as] In + 2A W{a - 2)1 + VI ] , 

a — A J 



= ^ !^-Ps^^Hn^ + 2Xiau;^ + s^^^)j, (49a) 
= 0, 

Substituting these expressions into the Baldin SR, Eq. (6), we obtain: 

eV fro/-, A 2 o 4N 1 21, 406 - 737/^2 + 304/x4 - 36//6 



{aE + /3M)r°) = I [3(1 - 4/.^ + 2/) + i/x^] In + 



6(4-^2)2 



44 - 788u2 + 1500u^ - 899ri^ + 215uS - 18ui° R ^ , , 

+ 3M4-, 2)5/2 ^-tg^- - 1 ^ (50a) 



(«. + /3M)r«)=^|^{ln/. 



1 

+ 



(4-/X 



2^2 



3/xV4^^ V 



(50b) 



Note that an identical result is obtained in the conventional one-loop Feynman diagram 
calculation [17, 18]. The corresponding (chiral) /i-expansion reads: 

n ^(RLO^ eV 11 f. 48(4 + 31n/x) 1521 o 1 .... ^ 

/ ^ N^iJLO^ eV 11 f, 4(1 + I21n/X) 117 , 1 rr■^^^ 

or, numerically (using g'^/Ai: = 13.8, M = 0.9383 GeV, ^l = 0.148), 

{oiE + PmYp^^^^ = 14.5 - 5.2 - 5.5 + ... = 3.8, (51c) 

(ae + /^m)!'^^^^ = 14.5 - 5.4 - 0.4 + ... = 8.7, (51d) 

in units of 10~"^ fm^. As one can clearly see, the fully relativistic leading order result is 
substantially different from the non- relativistic (heavy-baryon) limit. The higher order 
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corrections, which are suppressed by m^r/Af ~ 1/7, and hence are expected to be small, 
appear with large coefficients and generate a substantial modification of the leading order 
result. This (relativistic) efi^ect is likely to allow one to accommodate the relatively large 
A-resonance contribution to the magnetic polarizability [19]. 

For the forward spin polarizability, Bernard, Hemmert and MeiBner obtain[20] : 



7o = 



2171II 



IGtt^/x^ 

while, using the sum rule of Eq. 8, we obtain : 



(20 + 26 In fj,) f/ + 



1875 TT/x^ 
64 



+ ... 



(52) 



7o 



167r^/x^ 



1 - 



21 71 fl 

8 



59 



+ 26 In /X + 



1875 71 fi^ 
64 



(53) 



The difference between both expressions is confined to the analytic terms. The sum rule 
calculation has analyticity built in explicitly. 



VII. CHIRAL EXTRAPOLATIONS 

It is instructive to examine the chiral behavior of the one-loop result for the nucleon 
magnetic moment. Expanding Eq. (43) around the chiral limit (m,r = 0) we have 

= ^|l-27r/.-2(2 + 51n/.)/.2+^^3 + (5(^4)|^ (54^) 

= ^|-4 + 27r/.-2(l-21n/.)/.^-^/.=^ + 0(/.^)}. (54b) 

Apart from the first term, the constant, which renormalizes the counter term as described 
above, this expansion corresponds with the heavy-baryon expansion. The term linear in 
pion mass (recall that n = mT^/M) is the well-known leading nonanalytic (LNA) correction. 
On the other hand, expanding the same expressions around the large limit we find 



4'°°'^ = yfr^ (5 - 4 In /.) ^ + 0(//-^) , (55a) 



> - (|^(5-41n;.)-l 
''""'^ = 2(3 - 41n/x)-i + 0{ix-'). (55b) 

What is intriguing here is that the one-loop correction to the nucleon a.m.m. for heavy 
quarks behaves as l/rriquark (where rriquark ^ ^V)-! precisely as expected from a constituent 
quark- model picture [21]. Here this is a result of subtle cancellations in Eq. (43) taking place 
for large values of m.,^. In contrast, the infrared regularization procedure [14] gives the result 
which exhibits pathological behavior with increasing pion mass and diverges for = 2M. 

Since the expressions in Eq. (43) both have the correct large and small behavior they 
should be better suited for the chiral extrapolations of lattice results than the usual heavy- 
baryon expansions or the "infrared-regularized" relativistic theory. This point is clearly 
demonstrated by Fig. 5, where we plot the m7r-dependence of the full [Eq. (43)], heavy- 
baryon, and infrared-regularization[14] leading order result for the magnetic moment of the 
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proton and the neutron, in comparison to recent lattice data[22]. In presenting these results 
we have added a constant shift (counter-terms ko) to the magnetic moments, i.e., 



I^p = (l + «op + 4'°°P))(e/2M), 
l^n = («:on + 4'°°'^)(e/2M) 



(56) 
(57) 



and fitted them to the known experimental value of the magnetic moments at the physical 
pion mass, /ip ~ 2.793 and /x„ ~ —1.913, shown by the open diamonds in the figure. For the 
value of the nNN couphng constant we have used g'^/^n — 13.5. The m^-dependence away 
from the physical point is then a prediction of the theory. The figure clearly shows that 
the SR-motivatcd extrapolation, shown by the dotted lines, is in a good agreement with the 
mTT" dependence obtained in lattice gauge simulations. 

It is therefore tempting to use the SR results for the parametrization of lattice data. For 
example, we consider the following elementary two-parameter form: 



where kop and kon are fixed to reproduce the experimental magnetic moments at the physical 
m7r. The parameter a can be fitted to lattice data. The sohd curves in Fig. 5 represent the 
result of such a single parameter fit to the lattice data of Ref. [22] for the proton and neutron 
respectively, where = 1.6/M^ and a„ = 1.05/M^, and M is the physical nuclcon mass. 
This parametrization incorporates the experimental value at the physical pion mass and 
reproduces very well the trend of the lattice data. 

We would like to stress here that at larger pion masses the ChPT result should be 
interpreted as merely a convenient parametrization of lattice data, not as a prediction of the 
theory. In this way one obtains a parametrization which fits the lattice data at larger values 
of and has the correct dependence at lower values of m-;^. Apparently, the relativistic 
ChPT result which is consistent with analyticity provides the most convenient ground for 
the fit to lattice data, because it has an l/rUg behavior at larger values of m-^. It is useful to 
test the consistency of lattice simulations at larger and the experimental values based 
on such a single-free-parameter form which encodes the correct chiral behavior at low pion 
masses. 

The importance of the relativistic effects is even more emphasized in the polarizabilities. 
In Fig. 6 wc have plotted the ratio of one- loop Qe + Pm, Eq. (50), to its leading non- 
analytic term [first term in the expansion Eq. (51)]. Recently, the electric and magnetic 
polarizabilities of hadrons have been calculated in lattice QCD using the external field 
method [23, 24]. This method amounts to extract the electric and magnetic polarizabihties 
of hadrons from the quadratic term in the mass shift in an external electromagnetic field. 
The low statistics, pioneering results have so far been obtained for pion masses above 450 
MeV. As higher statistics calculations will become available in the near future [25, 26], our 
sum rule calculation of the pion mass dependence of the forward polarizabilities provides a 
way to connect these next generation lattice results to experiment. 




(58b) 



(58a) 
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VIII. PROSPECTS FOR THE FUTURE 



In this section we outline a procedure for two-loop calculations of the a.m.m. by using 
derivative forms of the GDH SR. Again, the GDH SR in its original form is practically 
unsuitable for such calculations, because in order to compute k in, e.g., QED to order 
one needs to compute the photoabsorption cross-sections to order — also to two loops! 
By contrast, using the derivative trick of Sect. Ill, we will need the (derivatives of) cross- 
sections up to order ctg^ only — one loop forms. 

Introducing the 'trial' value of kq, such that the full a.m.m. is defined as k — kq + Sk, 
and the GDH cross-section can be written as 



A(7(ko) = Act Act' + Aa" + ... , 
the first and second derivatives of the GDH SR at kq — read, respectively: 



5k + 5k 5k' 



M2 



— Act' 



25k' + {5k'Y + 5k5k" = 



I 



duj 



UJ 



Aa" 



(59) 

(60a) 
(60b) 



»;2 • 

em 



Now, consider Eq. (60a) at order al 

^TT^aem J OJ 



(61) 



where the superscript indicates the order of to which the quantity is considered. We 
can use Eq. (60a) and Eq. (60b) at order to obtain, respectively: 



5k^^^ = 



5k'^^^ - 



47r2 



oo 

J ^ 



duo 



UJ 



(62a) 
(62b) 



Substituting these in Eq. (61) we find the linearized GDH sum rule at order a 



2 . 
em" 



M2 



^Aa'(^)-^ 

UJ Sn'^ap, 



' duj 



UJ 



Aa'(2) 




Aa"(2) 



(63) 



and, using this sum rule one can obtain all two-loop corrections to the a.m.m. by computing 
the first derivative of the GDH cross-section Aa' at one-loop and tree level, and a second 
derivative of the GDH cross-section Aa" at tree level. 

For the QED case we have thus far worked out the required derivatives at tree level only. 
The first derivative is given in Eq. (26), while the second derivative is given by: 



A<7"(^) 



9 + 29X+ 18^2 - lOx^ / 9 , , , 

+ I 1 + — 1 log(l + 2x) 



2x 



(64) 
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with X — LojM. As already shown, the dispersion integral over the first derivative yields the 
Schwinger correction, c/. Eq. (27). The integral over the second derivative is divergent and 
to evaluate it we introduce an ultraviolet cutoff A 3> M: 



The same cutoff needs to be introduced in the first integral of the sum rule Eq. (63), such 
that the total result is independent of the cutoff. Thus, in QED we obtain 

5.(^) = -4^ /^A.'(3)-i%f51og^-6) . (66) 



In order to complete this two-loop calculation of the electron a.m.m. in QED one would 
need to evaluate the first derivative Act' at order which is given by the interference 
of the tree-level and one-loop Compton scattering ampfitudes, as well as by the tree-level 
bremsstrahlung and pair production mechanisms, all with one insertion of the Pauli vertex. 
Diagrammatically this calculation is depicted in Fig. 7. It is a future challenge to perform 
such calculation and to verify whether the linearized GDH SR reproduces the two-loop result 
obtained by usual techniques [27]: 

^ -0.328479 {alJi^) . (67) 

This would be an exact test of the GDH SR at the two-loop level. 

In the case of the nucleon, sum rules analogous to Eq. (61) may provide a more effi- 
cient method to do a two-loop calculation for the nucleon a.m.m.or for the forward spin 
polarizability, since much is already known about the one-loop ChPT amplitudes of pion- 
photoproduction [28], which are required in such a calculation. 



a: 



197 
144 



TT 



+ y(i-log2) + ^C(3) 



IX. CONCLUSION 

Direct application of QCD to low energy hadronic physics is made challenging by the 
feature that the quark/gluon degrees of freedom in terms of which QCD is couched are 
confined within hadronic systems. Nevertheless, the chiral EFT of QCD — chiral pertur- 
bation theory — allows a predictive description of the low-energy hadronic reactions. Also, 
by focusing on the chiral structure of the underlying Lagrangian such EFT methods per- 
mit one to make a link to lattice QCD calculations, even in the case where the mass of 
the underlying Goldstone bosons is considerably heavier than found experimentally. These 
are the two fronts which at present make the chiral EFT indispensable in relating QCD 
to low-energy observables. On the other hand there remain significant problems. Indeed, 
most such calculations involving nucleons are carried out at low orders in so-called heavy 
baryon chiral perturbation theory, which involves an expansion in powers of the inverse nu- 
cleon mass. On the other hand, we have demonstrated above by simple examples involving 
the nucleon magnetic moment and polarizabilities, that manifestly relativistic calculations, 
which include chiral corrections to all orders in = rriT^/M]^, do a better job than the 
"heavy-baryon" ones on both fronts in capturing the full implications of the theory. 
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For this reason it is important to have at one's disposal methods which allow such all 
orders evaluation of experimental quantities. This can in principle be achieved using conven- 
tional relativistic chiral loop calculations, but the price in terms of the number of diagrams 
which must be evaluated is high. For example, a one-loop evaluation of the nucleon polar- 
izability involves fifty- two such diagrams for the proton and twenty-two for the neutron [29]! 
Above we have presented an alternative approach, involving the use of dispersion relations 
and their related sum rules. The calculations done in this work were based on real-Compton- 
scattering sum rules, such as GDH and Baldin sum rules. However, the results are identical 
to what one would obtain in the usual loop calculations, provided no manipulations, such 
as infrared regularization which change the analytic structure are made. As is shown in the 
works of Gegelia et al. [13], there is no problem with power counting in this straightforward 
formulation of covariant ChPT, if the renormalization is performed in a suitable way. Both 
the infrared-regularized and the straightforward formulation of relativistic baryon ChPT 
resum the nominally higher-order relativistic effects, albeit differently. We stress that the 
resummation in the straightforward version is done in accordance with the principle of an- 
alyticity. 

We demonstrated above the utility of taking derivatives of the GDH sum rule, in order 
to convert it to forms which are sometimes more calculationally robust. We showed that in 

some cases these modified versions of the sum rules are equivalent to sideways dispersion 
relations, which involve taking the final (or initial) nucleon four-momentum off-shell. Using 
such relations, we straightforwardly evaluated fully relativistic nucleon magnetic moments 
and polarizabilities, using only tree level inputs and showed that they were fully equivalent 
to values obtained via much more labor-intensive loop techniques. We showed that such 
methods should also permit a future evaluation of the two loop electron anomalous magnetic 
moment. 

The successes demonstrated above suggest an obvious direction for future work. In par- 
ticular, the use of dispersive methods would seem to be ideal for a simplified way to perform 

two-loop evaluations. An obvious first test would be to use one-loop corrected photoproduc- 
tion amplitudes in order to produce corresponding one-loop photoabsorption cross sections, 
which can in turn be used to generate two-loop values of polarizabilities and of magnetic 
moments which are relativistically correct. These can then be compared to ongoing calcula- 
tions at two loop in the heavy baryon expansion for the forward spin polarizability and with 
the "experimental" value obtained from the sum rule. In conclusion there is much which 
can be done with such methods, and many future applications can be envisioned. 
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APPENDIX A: CALCULATION OF PHOTOPRODUCTION CROSS-SECTIONS 



We write the pion photoproduction amplitude in the following general form: 



7;^ = u{p') J2 Ms, t) Oi{e) u{p) , 



(Al) 



i=l 



where A's are the invariant amplitudes and O's are given by 

01 = 75, 

02 = £-7-?75, 

03 = 7 -£75, 

04 = 7-?75, 

with£-7-g = ^['j^^ju] s^q"- 

The total and double-polarized cross-sections are given respectively by 



(A2) 



Act 



27rA 



X 



27rA 



X 



4 /•! 



f d{cos9)YijAiAj, 

i,j=l -^-^ 
4 .1 

/ (i(cos 9) Ai Aj, 



(A3) 

(A4) 



with the scattering angle and 

y,, = IV[A(+)(/)0,A(+)(^5)0,] , 

A(+)(/)0,AW(j^)E-p0J , 



(A5) 
(A6) 



where S = 70757, p is the 3-momentum of the incoming nucleon, and the projection operator 
A(+^ is defined in Eq. (30). More explicitly, the traces are given by 



Y 



( \t -xi i y-x\ 

—x^ —2xy X + y 
^ x + y -2 + \i -xi 
— X —x^ —2xy j 

/ xi -e \ 

x^ 2xy —X — y 

-i -x-y 2 - it 



V 







(A7) 



(A8) 



xi j 



where x = {s - M^)/2M^ = u/M, y = (M^ - u)/2M^, ^ = p' ■ e = A^^ sin^, s = s/M"^, 
t = t/M'^ (s, t and u are the Mandelstam variables). 

In the notations of Eq. (37) the expression for the total cross-section of pion photopro- 
duction in Born approximation, including the nucleon a.m.m. terms, takes the following 
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form: 



.(TrOp) 



(TT+n) 



+ 2A [x'^{a - 2)(1 + Kpf + s/i'^ + ^nlas] } , 



27rC 



(3 s jj? In 



(3 + \ 



t; ..2\ 



+ 2A (aa; + s n 



I X 



r(.-^ P) 



ttC 
2nC 



X In 



a 


+ A 


a 


-A 


a 


+ A 


a 


-A 


a 


+ A 


a 


-A 



X{2 - as + 2{1 + x)p - f/) 



+ [-4 + (2 + s)q;]A|, 

P + \ a + \ 



H\sp-ii') In 



/5 — A a — A 



+ 2sf/X 



+ X^{Kp + Kn) {2X - In ^3^^ + + - ^f/) 



I X ^77. 



A,\ ^ 1 , a + A 
;i + /3)A+-ln— — 
z q; — A 



(A9) 



APPENDIX B: FIRST DERIVATIVES OF GDH BORN CROSS-SECTIONS OF 
PION PHOTOPRODUCTION 



In this appendix we present explicit expressions for the Aui cross-sections for the pion 
photoproduction on the proton. 



1. Pseudo- vector coupling 

Computing the Born graph contribution with the pseudo- vector ttNN couphng we obtain: 

Aai;°^) = .^l(2as+l-x)\n^^-2X[x{a-2-i^^) + s{2 + a)]\ 
^ M^x^ 1^ a — X J 

= A.(.»^) + |gA,= (Bl) 
A^<;°" = = = (B2) 

= l£ f^I + - ^«)} . (B3) 

A.ir»l = A<")-|^c.-A. (B4) 

where the notations are detailed in Eq. (37), and Acr*^'^"^^ is given in Eq. (36). 
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2. Pseudo-scalar coupling 

In the case of pseudo-scalar ttA^A^ coupling the expressions are different. 

= {2A 3 + ^^') + (3 - ^^^) In , (B5) 



nC ( 


2A {a- 


M^x^ { 










~ M^x' 



Aa\: = AaZn" = ^^Zn' = (B6) 
Aaif") = Aa[l^-^ = ^XP. (B7) 

Note, however, the contribution of these cross-sections to the linearized GDH sum rule 
Eq. (41) is exactly the same in both the pseudo-vector and pseudo-scalar cases. 

3. Integrands of the sideways dispersion relation 

The integrands entering the sideways dispersion relations of Eq. (44) for the proton a.m.m. 
are given by : 

A^S;-' = A^lf "', (B8) 

APPENDIX C: SECOND DERIVATIVES OF GDH BORN CROSS-SECTIONS 
OF PION PHOTOPRODUCTION 

In this appendix we list the expressions for the second derivatives of the GDH cross 
section, defined as: 



A ^ ^ A / 

Using the pseudo-vector coupling, we obtain: 



K0p=K0n=0 



(CI) 



= = ^ I A(4 - 2x/3 - 2a + /.^) - In -— ^ ^ , (C2) 

= a4;1 = ^ I A[2(.- + x)a -2x + ^i^]-ln^\, (C3) 
M _ A^W _ /o\r oaf, , o^„2^ - " + A 



Aa!fJ^ = Aa'^ = ^ | 2A[-2/5(l + x) + - 2 + /.^ + In -— ^ } . (C4) 



Using the pseudo-scalar coupling, we obtain: 

lA.<-'> = A<"' = -^{2A(l+;5)-ln^ 

= (m' - 2./J) = Ha (/J - ax) , (C6) 
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FIG. 1: Tree-level Compton scattering amplitude in QED. 




FIG. 2: Tree-level Compton scattering graphs with one Pauli-coupling insertion. The circled vertex 
corresponds to the Pauli coupling. 




FIG. 3: Tree-level pion photoproduction graphs. The circled vertex corresponds to the Pauli 
coupling. 
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FIG. 4: Absorptive part of the jNN vertex with final nucleon off the mass shell, used in the 
sideways dispersion relation for the nucleon a.m.m.. Diagram (a) : TT~^n loop where the photon 
couples to the vr"'"; diagram (b) : n^p loop where the photon couples to the charge of the proton. 
The vertical dotted lines indicate that the vrA^ intermediate state is taken on shell. 
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FIG. 5: Chiral behavior of proton and neutron magnetic moments (in nucleon magnetons) to one 
loop compared with lattice data. "SR" (dotted lines): our one- loop relativistic result, "IR" (blue 
long-dashed lines): infrared-regularized relativistic result, "HB" (green dashed lines): LNA term 
in the heavy-baryon expansion. Red solid lines: single-parameter fit based on our SR result. Data 
points are results of lattice simulations. The open diamonds represent the experimental values at 
the physical pion mass. 
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FIG. 6: The ratio of the one- loop proton (solid curves) and neutron (dashed curves) forward 
polarizabilities to their LNA terms in the heavy-baryon expansion. 
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FIG. 7: Interference contributions which give rise to the Aa'^^^ in QED. The circled vertex corre- 
sponds to the Pauh couphng. 
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